Nonlocal Nonlinear Optics in cold Rydberg Gases 



S. Sevincli, 1 N. Henkel, 1 C. Ates, 2 and T. Pohl 1 

1 Max Planck Institute for the Physics of Complex Systems, 01187 Dresden, Germany 
2 School of Physics and Astronomy, University of Nottingham, Nottingham, NG7 2RD, United Kingdom 

We present an analytical theory for the nonlinear optical response of a strongly interacting Ry- 
dberg gas under conditions of electromagnetically induced transparency. Simple formulae for the 
third order optical susceptibility are derived and shown to be in excellent agreement with recent 
experiments. The obtained expressions reveal strong nonlinearities, which in addition are of highly 
nonlocal character. This property together with enormous strength of the Rydberg-induced non- 
linearities is shown to yield a unique laboratory platform for nonlinear wave phenomena, such as 
collapse-arrested modulational instabilities in a self-defocussing medium. 
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Advances in designing materials with highly intensity- 
dependent refraction [1-3] have ushered in numerous 
studies of nonlocal nonlinear wave phenomena [4-8]. 
Many of these settings, such as, nematic liquid crystals 
[1, 2] or thermal media [3], require high power laser light. 
On the other hand, electromagnetically induced trans- 
parency (EIT) in ultracold multi- level atoms [9, 10] pro- 
vides an elegant mechanism to suppress photon loss and 
simultaneously increase light-matter interaction times to 
enhance nonlinear effects. Combined with sufficiently 
large nonlinearities, this holds great potential for few- 
photon nonlinear optics [11, 12] and may enable appli- 
cations in communication and quantum information sci- 
ence. 

Recently, it was recognized that EIT-schemes involving 
highly excited atomic Rydberg levels provide promising 
perspectives for such applications [13-22]. In particular, 
the huge polarizability of Rydberg states gives rise to 
giant Kerr coefficients [15], but also entails strong long- 
range interactions, which render Rydberg-EIT media in- 
trinsically nonlinear. Indeed, a recent theory for two- 
photon pulses revealed the emergence of strong effective 
photon-photon interactions [23], while experiments [21] 
and numerical calculations [24] demonstrated greatly en- 
hanced nonlinear absorption coefficients in the opposite 
limit of large photon numbers. 

In this letter, we develop an analytical theory for 
the nonlinear optical response of a strongly interacting 
Rydberg-EIT medium to monochromatic multi-photon 
light sources. Based on the approach, we give a simple 
formula for the nonlinear absorption coefficient that pro- 
vides an excellent description of recent measurements on 
cold Rubidium gases [21]. For large single-photon detun- 
ings, absorption is shown to be greatly suppressed - yet 
maintaining huge refractive nonlinearities, that exceed 
previous records in ultracold Kerr media [10] by several 
orders of magnitude. Combined with their long range 
this makes for an ideal nonlinear medium to study non- 
local wave phenomena, in which the strength, the range 
and even the sign of the nonlocal interaction kernel can 
be widely tuned with high accuracy. To demonstrate this 




FIG. 1: a) Three-level scheme for isolated atoms, where the 
atomic ground state |1), an intermediate state |2) and a highly 
excited Rydberg state |3) are mutually driven by a strong con- 
trol and a weak probe field with Rabi frequencies Q c and Q p , 
respectively. On two-photon resonance, EIT ensures lossless 
propagation of the optical fields, unaffected by spontaneous 
decay (7) and the single-photon detuning A/2, (b) In a gas 
of atoms, the strong van der Waals interaction between atoms 
in Rydberg states (|3)) inhibit multiple Rydberg excitations 
within a blockade radius R c , giving rise to a strongly non- 
linear optical response of the medium. The resulting non- 
linear beam propagation, for example, leads to modulation 
instabilities, as shown in (c) for a Rubidium 705*1/2 Rydberg 
gas with a density of 8 x 10 13 cm~ 3 and Q p /2tt = 0.35MHz, 
Q c /2tt = 80MHz, A/2tt = 1.2GHz. 

potential, we present numerical results for the propaga- 
tion of cw laser light and show that paradigm phenom- 
ena, such as optical solitons and modulational instabili- 
ties (see Fig.lc) are observable with current experimental 
capabilities. 

Consider first the propagation of a beam with ampli- 
tude Q p (see Fig. lb) and wavenumber k as described by 
the paraxial wave equation 

where Vi accounts for the transverse dynamics with re- 
spect to the axial coordinate = (x, y) perpendicular 
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to the propagation direction z. The relevant medium 
properties are contained in the complex susceptibility 
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which is determined by the dipole matrix element pi2 of 
the probe transition and the corresponding atomic coher- 
ence density pi2- The probe field, fi p , drives the lower 
transition between the ground state, |1), and a low-lying 
excited state, |2), of ladder- type three- level atoms (see 
Fig. la), whose optical response is controlled by a strong 
control field, driving the upper transition between |2) and 
a Rydberg state, |3), with a Rabi frequency Q c > Q p . 
Without interactions, this yields a perfect EIT medium, 
in which each of the N atoms in the gas settles into a 
dark state \di) oc Q c | 1^) — Q p (r^)|3i) (i = l,...,iV) such 
that pi2 = X = an d the probe beam is unaffected by the 
atomic medium [9]. In the presence of strong Rydberg- 
Rydberg atom interactions the gas dynamics becomes 
highly correlated due to the resulting level shifts of mul- 
tiply excited Rydberg states. Within a critical blockade 
radius R c all but a single Rydberg excitation are inhib- 
ited [25] (see Fig. lb) and removed from two-photon reso- 
nance, thereby diminishing EIT, and, thus giving rise to 
nonlocal absorption and refraction within a range ~ R c . 
Since the Rydberg state population in the unperturbed 
dark states \di) is proportional to ^ p (r^) 2 one, hence, 
expects an intensity-dependent, i.e. nonlinear, optical 
response. 

Having established a simple picture of the basic mech- 
anisms we now derive the resulting optical susceptibility 
from the underlying Heisenberg equations for the atomic 
transition operators <r^ = |a^)(ft| (a,/3 = 1,2,3). In 
the limit of low probe intensities (f2 p (r^) <C fl c ) these 
can be expanded in Q p /Q, c [26]. Upon adiabatic elimina- 
tion of one obtains a single dynamical equation for 
the two-photon transition operator of the ith atom 



dt* 13 ~ 



Q p (rj) + fl c <3i3 7i3 



where 033 = 03^013 , T = 7 + 712 — iA, A/2 is the single- 
photon detuning and the rates 7, 712 and 713 account 
for the spontaneous decay of the intermediate state as 
well as the linewidth of the probe and two-photon tran- 
sition, respectively. The last term in eq.(3) describes the 
interactions between atoms in the Rydberg state |3) and 
Vij — Cq 1 1 — Tj\ 6 denotes the corresponding van der 
Waals potential for atoms at positions and rj. Since 
the van der Waals coefficient Cq oc n 11 drastically in- 
creases with the atom's principal quantum number n, 
the interaction between highly excited Rydberg atoms 
exceeds that of the two low- lying states by many orders 
of magnitude. Proper inclusion of the resulting strong 
atomic correlations requires knowledge of the two-body 
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FIG. 2: Nonlinear transmission of a cold Rubidium Rydberg- 
EIT medium with |3) = \60Si/ 2 ) at two different densities 
and control Rabi frequencies, and for 712/2^ = 110kHz and 
7i3 /2tt = 220kHz [21]. Up to Q p0 « 0.3Q C there is good 
agreement between our low-Q p prediction eqs.(5) and (6) 
(solid line) and the experimental data [21, 27] (symbols). The 
dashed lines neglect the drop in absorption due to attenuation 
and averaging over the initial transverse beam profile. 



correlators o^o~^^ whose dynamics follows from eq.(3) 
by applying the chain rule. Being primarily interested 
in the leading order nonlinear contribution to x> we can 
once more expand the resulting two-body equations to 
leading order in Q p . This amounts to dropping direct 
three-body correlators and, thus, yields a closed set of 
evolution equations for the one- and two-body opera- 
tors. Setting ±{a%) = = 0, the steady 

state expectation values are then readily obtained from 
the resulting set of algebraic equations. Finally, we 
take the continuum limit by defining continuous densi- 
ties p a p(r) = Ei(^)^( r - r i)> and obtain 
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where p = S(r — r^) is the total atomic density. To- 
gether with eq.(2) this yields the leading-order nonlinear 
susceptibility and permits to propagate the probe beam 
according to eq.(l). 

If the atoms are driven on single-photon resonance, 
A = 0, the main interaction effect will be nonlinear ab- 
sorption. Hence, one can neglect the transverse beam 
dynamics (Vj_) as well as the nonlocality in eq.(4), by 
setting ^p(r') « ^p(r). With this simplification one ob- 
tains local first and third order susceptibilities, defined by 
x(r) = + X^ 3 ^p( r )- The remaining spatial integral 
in eq.(4) can be carried out analytically to give 
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FIG. 3: (a) Effective photon-photon interaction potentials, 
introduced in eq.(8). (b) Growth rate Tmi of intensity mod- 
ulations with wavenumber k for defocussing nonlinearities of 
different strengths Q 2 a. The dashed lines show the corre- 
sponding imaginary part while the solid lines correspond to 
the real part of Tmi- The critical value of Q 2 «mi = 50.06 
marks the onset of a modulational instability within a narrow 
window of wave numbers indicated by the grey shaded area 
for Q 2 a = 80. 



This expression permits a simple interpretation, by in- 
troducing the resonant blockade radius R C1 defined by 
the distance at which the interaction \Cq\/R® exceeds 
the width Seit = of the EIT window [23]. Sub- 

stitution of C 6 by R c = (|C 6 |/(5eit) 1/6 shows that x (3) 
is proportional to the corresponding two-level response 
times the number R 3 p of blockaded atoms, which is con- 
sistent with the simple picture outlined above and the 
numerical findings of [24]. 

Experimentally, nonlinear absorption has been re- 
cently studied in a cold Rubidium gas involving |3) = 
|605i/ 2 ) Rydberg states [21]. In the experiments the 

transmission, T, of a Gaussian probe beam (Q p = 

— 2 / 2 \ 

Q p oe~ r±/w ) through the gas of length I was measured 
for different intensities and atomic densities. Within the 
local approximation this configuration permits a simple 
solution of eq.(l) for the integrated beam transmission 
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where p = ft£ oXl (3) (l - T )/xi (1) and T = e~ k ^ 1)l is 
the first order transmission. Fig. 2 shows a comparison 
to the measured transmission for two different densities 
and demonstrates good agreement, even for rather large 
probe Rabi frequencies of up to £7 p o ~ 0.3^ c . Note that 
the backaction of the nonlinear beam attenuation onto 
susceptibility and, equally important, the averaging over 
the transverse beam profile are both essential for a proper 
description of the experiment. Neglecting these effects 
yields the dashed lines in Fig. 2, which significantly over- 
estimates the nonlinear absorption. 

Since on resonance xj 3 ^ = IXr^I e Q-(5)], large non- 
linear refraction is inevitably accompanied by high pho- 
ton loss. However, for large single-photon detunings A ^> 
7 eq.(4) yields xi °° (t/A)Xr , such that dissipative loss 



can be greatly suppressed. For instance, for a Rubidium 
Rydberg gas with ^ c /2tt = 5MHz, ^ p /2tt = 0.5MHz, 
A = 30GHz and p — 8 x 10 13 cm -3 one obtains a large 
absorption length of / a b s ~ 1mm, and yet a high nonlin- 
ear refractive index ri2 ~ 2 x 10 4 cm 2 /W which is 5 orders 
of magnitude greater than previously obtained with ul- 
tracold Rb groundstate atoms at the same density [10]. 

As refraction starts to dominate absorption, the nonlo- 
cality of x^ [cf. eq.(4)] becomes significant. To account 
for its effects on the transverse beam propagation we re- 
cast eqs.(2,4) into 
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where we assumed 712 <C 7, 713 <C 5eit = ^c/A 
and introduced the off-resonant blockade radius R c = 
(C 6 /^eit) 1/6 (C 6 A > 0) [23] set by the off-resonant 
EIT width Seit- To simplify matters, we proceed by 
defining scaled coordinates r = z/(kR 2 ) : £ = r±/R c 
and the dimensionless probe amplitude ^, normalized to 
J f2 2 (£,r)<i 2 £ = 1. Retaining the local approximation 
along the propagation direction [31] this yields a two- 
dimensional nonlinear Schrodinger equation 
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the nonlinearity, P p denotes the probe beam power 
and the effective interaction potentials J7 m (£) = 
dz [l + (£ 2 + z 2 ) 3 ] m are shown in Fig. 3a. As a ex 
Cq, repulsive atomic interactions lead to self-defocussing 
nonlinearities, while attractive atomic interactions map 
onto self-focussing nonlinearities. 

The former case, can, e.g., be realized with cold 
Rb(n5 , 1 / 2 ) Rydberg states as in the experiments [21] 
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FIG. 4: Output beam profile |^| 2 (r ± ) for Q p0 /2tt = 0.35MHz, 
Q c /2tt — 80MHz and A/2tv — 1.2GHz after propagation over 
I = 210/im through a Rb(705i/ 2 ) EIT medium at three dif- 
ferent densities (a) 4 x 10 13 cm~ 3 , (b) 5.5 x 10 13 cm~ 3 and (c) 
8 x 10 13 cm -3 . For the color coding each distribution has been 
normalized by the actual maximum intensity |n max | 2 - 
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FIG. 5: Intensity profiles |^| 2 (r ± ) for a Sr(50 1 5 f ) EIT 
medium with Q /2tt = 0.3MHz, Q c /2tt = 15MHz, A/2tt = 
3.2GHz at two different densities of (a,c,e) 8 x 10 11 cm -3 and 
(b,d,f) 1.2 x 10 12 cm -3 . Panels (a,b) show the stable soliton 
solutions, while panels (c-f) show the compressed output in- 
tensity profile after a propagation length / = 240/xm for an 
input beam with (c,d) v — 2 and (e,f) v — 6. The color coding 
is identical to Fig. 4. 



discussed above. While the resulting photon-photon in- 
teractions are isotropically repulsive, the corresponding 
momentum space interaction U\{k) is not sign-definite. 
Despite being defocussing, the present nonlocal inter- 
action can, consequently, promote a modulational in- 
stability. In Fig. 3b we show the corresponding rate 

Tmi(&) = -|\/^ 2 - 4a^ 2 Z7i(&) [4] for a given mode with 
wave number k to grow out of a homogenous amplitude 
f2(£) = const.. This growth rate Tmi assumes real values 
within a narrow range around k ~ 2tt at a critical interac- 
tion strength aft 2 « 50 [28], resulting in stable transverse 
intensity modulations on a length scale ~ R c as the beam 
propagates through the medium (cf. Fig.lc). To examine 
their observability, we performed numerical simulations 
of eq.(8) for a Rb( 7051/2 ) Rydberg gas traversed by a 
super-Gaussian beam Q = ^ e -( ^ w ^ + ^ with v — 6 and 
small spatial phase noise <\>. Fig. 4 shows calculated out- 
put intensity profiles for different atomic densities and 
demonstrates that highly localized, rather regular inten- 
sity patterns can be realized in high density gases with 
feasible laser parameters. 

Self-focussing nonlinearities arise from attractive Ry- 
dberg interactions, as occurring between n 1 ^ states of 
Strontium atoms [29], for which EIT has been recently 
observed [14]. In this case, modulational instabilities can, 
in principle, occur for any a < and sufficiently large 
beam widths. More importantly, however, the attrac- 
tive photonic softcore-interaction enables the formation 
of stable bright solitons, leading to tight beam focussing. 
A simple variational analysis of eq.(8) (see, e.g., [30]) 
yields a critical interaction strength of a so ~ 0.71, above 
which stable bright solitons exist. Figs. 5a and 5b show 
two examples for a Sr(50 1 5'o) gas and reveal a character- 



istic soliton size < R c . In addition we show final inten- 
sity profiles for a Gaussian (y = 2) and super- Gaussian 
{y — 6) input beam. All cases are for an input width of 
w — 3R C and demonstrate significant focussing after the 
considered propagation length / = 240/im. In the latter 
case, the beam compression is superimposed by a radial 
MI leading to a ring-shaped hollow output beam. With 
a typical size of several /im these structures are readily 
observable experimentally. 

In conclusion, we have presented a theory for the non- 
linear response of a strongly interacting Rydberg-EIT 
gas, giving good agreement with recent measurements. 
The derived expressions for the third-order susceptibil- 
ity suggest that huge nonlinearities of highly nonlocal 
character can be experimentally realized, which provides 
an ideal setting to study complex nonlinear wave phe- 
nomena. To demonstrate these prospects we have shown 
that the observation of basic effects such as the forma- 
tion of bright solitons and collapse-arrested modulational 
instabilities are within experimental reach. The latter 
is particularly interesting in the uncommon case of de- 
focussing, nonlocal nonlinearities. Here, the repulsion 
between emerging intensity peaks combined with trans- 
verse beam confinement may promote the formation of 
transverse supersolid or crystalline states of photons. 
This question may be addressed within the present ap- 
proach, extended to quantum light in order to account 
for atom-photon and photon-photon correlations, which 
would open up a general framework for studying many- 
body physics with strongly interacting photons. From 
a different perspective, we expect the discussed nonlin- 
ear light propagation to be relevant for interpreting cold 
Rydberg gas experiments at high densities. 
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